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ABSTRACT 

To reconcile the observed unusual high luminosity of NuSTAR X-ray pulsations from M82X-2 
with the most extreme violation of the Eddington limit, and in view that the persistent X-ray 
radiation from M82X-2 almost precludes the possibility of common pulsars, we tackle the problem 
by the implications of microscopic theory of black hole (MTBH). The preceding developments of 
MTBH are proved to be quite fruitful for the physics of ultra-high energy (UHE) cosmic-rays. 
Namely, replacing a central singularity by the infrastructures inside event horizon, subject to 
certain rules, MTBH explains the origin of ZeV-neutrinos which are of vital interest for the source 
of UHE-particles. The M82X-2 is assumed to be a spinning intermediate mass black hole resided 
in final stage of growth. As a corollary, the thermal blackbody X-ray emission arisen due to the 
rotational kinetic energy of black hole escapes from event horizon through the vista to outside 
world that detected as ultraluminous X-ray pulsations. The M82X-2 indeed releases 99.6% of 
its pulsed radiative energy predominantly in the X-ray bandpass 0.3 — 30 keV. We derive a pulse 
prohle and give a quantitative account of energetics and orbital parameters of the semi-detached 
X-ray binary containing a primary accretor M82X-2 of inferred mass M ~ 138.5 — 226 Mq and 
secondary massive, M 2 > 48.3 —64.9 Mq, 0/B-type donor star with radius of i? > 22.1 — 25.7 Rq, 
respectively. We compute the torque added to M82X-2 per unit mass of accreted matter which 
yields the measured spin-up rate. 

Subject headings: black hole physics; accretion: accretion discs; X-rays: binaries; X-rays: individual 
(NuSTAR J095551-h6940.8) 


1. Introduction 

The physical nature of the accreting off-nuclear 
point sources in nearby galaxies, so-called the 
ultraluminous X-ray sources (ULXs), has been 
an enigma because of their high energy output 
characterized by super-Eddington luminosities up 
to two orders of magnitude higher than those 
observed from Galactic X-ray binaries, and un¬ 
usual soft X-ray spectra with blackbody emission 
around ;< 0.3 keV and a downturn above 5 


keV fe.g. Fabbiano et al. (19891): King 

et al. 

(I2OOIII: ISwartz et al. (2004h: iRoberts 

( 

2007 

); 

Gladstone et al (2009h: Feng & Soria 

( 

2011 

); 

Liul (201lt): Walton et al. 1 (2011): Liu et al. 1 

(2013): Motch et al 1 (2014)1. In spite of signifi- 


cant efforts in more than three decade since the 


discovery of ULXs, the astronomers have not yet 
clarified their nature. Without care of the histor¬ 
ical justice and authenticity, it should be empha¬ 
sized that the ULX sources have attracted a great 
deal of observational and theoretical attention, 
in part because their luminosities suggest that 
they may harbor intermediate mass black holes 
(IMBHs) with an ubiquitous featu re of the mass 
fits of more than 10^ — 10^ Mc^ (iFeng &: SorieT 


2011:[^lbert fc Mushotzkv Ill999l: Makishima et al. 


2000). A strong argument in favor of IMBHs is 
the presence o f a soft, 0.10.2 keV component in 


their spect ra ( Kaaret et al I 20031 : Miller et al. 


20031 120041) . The hyper luminous X-ray sources 
with luminosities > 10^^ ergs“^ are tho ught to be 


amongst the strongest IMBH candidates (jMatsumoto et al. 
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l2004l) . Assuming the emission is isotropic, in 
general, the extreme luminosities of ULXs sug- 


gest either the presence of IMBHs and sub- 
Eddington accretion (e.g. Komossa & Schulz 

( 1998h:lMiller & Colbert 

(20041: Dewangan et al. 

(2006ll: Stobbart et al. (2006 

1: Patruno et al. 

); 

(2006ll: Portegies et al. I( 

2 ()QiT); 

Kong et al. (2007 

Casella et al. (2008ll: 

Feng & Kaaret 1 (2010ll: 

Liu et a . (2013h: Fens' h Soria 1 (l201lh: Walton et al. 

( 2 OI 3 II: Pasham et al. (2014lll. or stellar-mass 


black holes (M < 10 Mq) that are either breaking 
or circumventing their Eddington limit v ia some¬ 


what g eometric beaming o f accretion flow (iBegelman 
2002: King et al. 2001: Okaiima et oi. I I 2 OOOI: 


Roberts 112007 : Zamoieri fc Roberts 20091: Stobbart et al. 


20061 : Poutanen et al. 2007ll . The latter remains 


still a relatively poorly understood regime. I Sutton et al. 
( 2012t l present the results from an X-ray and op¬ 
tical study of a new sample of eight extreme lu¬ 
minosity ULX candidates, which were selected as 
the brightest ULXs (with Lx > 5 x 10^°ergs“^) 
located within 100 Mpc identified in a cross cor¬ 
relation of the 2XMM-DR1 and RC3 catalogues. 

These objects seemed to be amongst the most 
plausible candidates to host larger, IMBHs. But, 
future improved observations are needed yet to de¬ 
cide on the issue. The quasi-periodic oscillations 
(QPOs) and frequency breaks in XMM-Newton 
power-density spectra of ULXs with luminosity. 

Lx > 10^°ergs“^, suggest that the black hole 
masses are more consistent with IMBHs than 


stellar-mass black holes (|Miller fc Colbert 1120041) . 


Or, using high sensitive dataset in the Fe K region 
obtained in scope of a long Suzaku program on 
Holmberg IX X-1 for any lumin ous, isolated ULX 
to date, IWalton et al. \ ( 2013 1 find no statisti¬ 
cally significant narrow atomic features in either 
emission or absorption across the 59 keV energy 
range. They conclude, therefore, that the mod¬ 
els of spherical super-Eddington accretion can be 
rejected, as can wind-dominated spectral models. 
The lack of iron emission implies that the stel¬ 
lar companion in sub-Eddington accretion onto 
an IMBH is unlikely to be launching a strong 
wind, and therefore the black hole must primarily 
accrete via roche-lobe overflow. There are also 
seve ral other new result s of m ore recent studies: 
e.g. Wiktorowicz et a/71 ( 20I5ll demonstrate that 
if the transferred mass is efficiently converted to 
X-ray luminosity (with disregard of the classical 


Eddington limit) or if the X-rays are focused into 
a narrow beam then binaries can form extreme 
ULXs with the X-ray luminosity of Lx ^ 10"*^ erg 
s“^. They envisage that these systems are not only 
limited to binaries with stellar-origin black hole 
accretors, but also can be identified with neutron 
star systems. For the latter the typical donors are 
evolved low-mass (2 Mq ) helium stars with Roche 
lobe overflow rate of ~ 10“^ MQyi~^. This study 
does not prove that any parti cular extreme ULX is 
a regular binary system. Or, iFragos et al. I (I 2 OI 5 II 
performed a population synthesis studies of ULXs 
with neutron star accretors. Combining binary 
population synthesis and detailed mass-transfer 
Jmodels, they conclude that the binary system that 
formed M82X-2 is most likely less than 50Myr 
old and contains a donor star which had an initial 
mass of approximately 8 — 10 Mq, while the NS’s 
progenitor star had an initial mass in the 8—25 Mq 
range. Or. iMushtukov et al. I ( 2015ll study prop¬ 
erties of luminous X-ray pulsars using a simplified 
model of the accretion column. The resulting lu¬ 
minosity of NS pulsars may can reach values of the 
order of erg s“^ for the magnetar-like mag¬ 
netic held, since the equilibrium, where the Alfven 
radius matches the co-rotation radius, indicates a 
magnetic held ol B > and long spin peri¬ 

ods {P > 1.5s). They conclude that a substantial 
part of ULXs are accreting neutron stars in binary 
systems. All these proposals have own advantages 
and difficulties. Nevertheless, no single theory has 
been invented yet which successfully addresses the 
solution to the problems involved. The aforemen¬ 
tioned studies are not exception to the rule that 
as phenomenological approaches they suffer from 
own difficulties. Namely, they are strongly model 
dependent, and subject to many uncertainties and 
controversies. The physics is obscured by multi¬ 
ple arbitrary assumptions and proliferation of a 
priori free parameters involved, while a consistent 
complete theory would not have so many free pa¬ 
rameters. So, the issue is still controversial, and 
these exotic results and conclusions are still lame 
and fo r sure not proven by the authors. For ex¬ 
ample, |wiktorowicz||el^iZ ( 2015 1 have evoked the 
possibility that the transferred mass is efficiently 
accreted onto a compact object and converted 
to X-ray luminosity in the full range of possible 
mass accretion rates, which is in contrast with 
the generally accepted view that the conversion 
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efficiency decre ases with increas i ng m ass accre¬ 


tion rate (e.g., iPontanen et al. ([20071)). Or, in 


case of model of the accretion column, a necessary 
condition for the X-ray luminosity to exceed the 
Eddington limit is a certain degree of asymme¬ 
try in t he distribution of matter over the Alfven 
surface ( Basko fc Sunvaev IQTtih . which is a qual¬ 
itative picture, perhaps, but not quantitative one, 
or a B > 10^^ G magnetic field of a magnetar- 
like pulsar is rather above the quantum limit, etc. 
We will not be concerned here with the actual 
discussion of these results, which require further 
extensive and careful analysis, in order one could 
verify the proposed estimates and constrain model 
parameters. There were still many open key ques¬ 
tions arisen inevitably that we have no under¬ 
standing of ULX physics. A main physical issue 
whether ULXs are powered by IMBHs or normal 
stellar black holes to date is unresolved, primarily 
because we do not have dynamical mass measure¬ 
ments of the compact objects that power ULXs. 
So, it is premature to draw conclusions and only 
time will tell whether any of these intriguing pro¬ 
posals is correct and which of the hypothesized 
ULX scenario is actually realized in nature. 

1.1. The ultraluminous X-ray pulsations 

The most striking is the recent revolutionary 
NuSTAR discovery of the first rare and mighty ul¬ 
traluminous X-ray pulsations with the maximum 
luminosity L(3—30keV) = 4.9xl0^®ergs“^, of av¬ 
erage p eriod 1.37 s with a 2.5 -daY sinusoidal mod¬ 
ulation ( Bachetti et al. llioill l. coming from an ul¬ 
traluminous source, NuSTAR J095551-I-6940.8, lo¬ 
cated nearby starburst galaxy M82 (NGC 3034). 
The pulsed emission centroid is spatially consis¬ 
tent with the location of a variable M82X-2 which 
further secures the association of the pulsating 
source with M82X-2. Detection of coherent pulsa¬ 
tions, a binary orbit, and spin-up behavior indica¬ 
tive of an accretion torque unambiguously, allow 
to feature the M82X-2 as mass-exchange binary 
that contains a nondegenerate secondary donor 
star. Eventually, in addition to the orbital modu¬ 
lation, about an evid ent linear spin-up of th e pul¬ 
sar was reported by (IBachetti et al. 11201411 . with 
p ~ —2 X 10“^°s/s over the interval from mod¬ 
ified Julian days 56696 to 56701 when the pulse 
detection is most significant. Phase connecting 
the observations enables detection of a changing 


secular spin-up rate over a longer timespan as 
well as erratic variations. Future observations will 
show whether the current spin- up rate is secu¬ 


lar. iKluzniak fc Lasota (120151 1 point out that 


the spin-up to luminosity ratio 10“®° (erg s)“^ is 
an order of magnitude lower than the typical ra¬ 
tio observed in the X-ray pulsars, which makes an 
interpretation of the data in terms of a strongly 
magnetized X-ray pulsar quite challenging. In ab¬ 
solute terms, this spin-up rate is orders of magni¬ 
tude higher than the values measu red in the usual 


accretion powered X-ra y pulsars (jBildsten et o.l. 
19971 : [ziolkowski 1985ll . 


1.2. Key objections 

A current understanding of NuSTAR discovery 
is quite complicated. At first glance it seems as 
though the NuSTAR team has demonstrated that 
the super-Eddington accretion is also possible in 
ULXs hosting a neutron star, because it is gener¬ 
ally believed that the pulsating X-ray sources are 
magnetic neutron stars whi ch are accreting matter 


from a binary companion (IPringle fc Rees 


Davidson fc Ostriker 1973r Lamb et al. 


ion (jP 
197.1 


1972: 


197.1: 


Lamb Il975f) . Therefore, it seems that there is 
nothing left but M82X-2, which until recently 
astronomers had thought was a black hole, is 
the brightest magnetic neutron star system ever 
recorded. This point of view is widely quoted 
in literature and, at first sight, seems eminently 
reasonable. However, deeper examination raises 
several disturbing questions, if the above result is 
really valid. Even though, if for a moment we take 
the classical accreting neutron star pulsar system 
as a basic assumption, note that this model is not 
flawless. The actual mechanism by which pul¬ 
sars convert the rotational energy of the neutron 
star into the observed pulses is poorly understood. 
Many theoretical models have been proposed that 
account for such features, but no sing l e one is 


198: 


comp elling (jSieber fc Wielebinski I Il98ll : iMichel 


Explaining periodic source M82X-2 that obvi¬ 
ously has black hole energetics with a ~ 1.4 Mq 
compact object using the accreting neutron star 
pulsar systems is extremely challenging, because 
of several problems we encountered. The difficulty 
becomes apparent when one follows the three ob¬ 
jections, which together constitute a whole against 
the claim that M82X-2, perhaps, is a common 
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pulsar: 1) The NuSTAR team discovery is the 
most extreme violation of the so-called Edding¬ 
ton limit, i.e. the pulsed luminosity of M82X-2 
reaches about ^ 26.9 times brighter than the the¬ 
oretical threshold at the spherical accretion for 
^ 1.4 Mq stellar-mass black holes where the out¬ 
ward pressure from radiation balances the inward 
pull of gravity of the pulsar. The accretion is 
inhibited once radiation force is equal or grater 
than gravity force. 2) The difficulty is brought 
into sharper focus by considering the associa¬ 
tion with M82X-2, which is featured with high 
luminosity (~ 1.8 x 10^°ergs“^) of additional 
persistent continuous broad X-r ay radiation ob¬ 


served earlier from its active state (|Bachetti et al. 


201411 . This more compelling argument in some¬ 
how or other implies the luminosity ^ 100 times 
if compared to the Eddington limit. Such a col- 
limation 100), which is comp arable to that 


obtain e d for black holes (e.g . . lOhsuea et al. 


(|2007ll : ISadowski fc Naravanf ( 2015ri l. would be 


needed to explain M82X-2 as beamed radiation 
from neutron star. 3) Equally noteworthy that 
the centroid of the persistent X-ray emission is be¬ 
tween M82X-2 and M82X-1. If M82X-1 is indeed 


harbors plausible IMBH, (e.g. 

Patruno et al. 

(20061; Portegies et al. (l2007^: Roberts ( 

2007 

); 

Casella et al. 

(I 2 OO 8 I; Feng & Kaaret ( 

2010 

); 

Feng & Soria 1 ( 

20111: IPasham et al. 1 (201^ 

1), we 


expect the similarity of the persistent X-ray prop¬ 
erties of the M82X-1 and M82X-2 to imply that 
the non-pulsed emission from the latter originates 
in the accretion disc, as it must in the black hole 
M82X-1. In this sense, the NuSTAR discovery is 
unexpected and still hard to be explained in the 
context of magnetic neutron star pulsar model. 
The fraction of ULXs powered by neutron stars 
must be considered highly uncertain and many 
details of this scenario remain poorly understood. 
Added to this was the fact that NuSTAR finding 
may indeed not be rare in the ULX population. 
In the future astronomers also will look at more 
ULXs, and it is a matter of time until they could 
prove an expected ubiquitous feature of even more 
energetic ULX pulsations, as being common phe¬ 
nomena in the Universe. If confirmed, this would 
support a scenario in which the more ULXs beat 
with the pulse of black holes rather than mag¬ 
netized neutron star systems. What if NuSTAR 
detection might radically change ones view of Na¬ 


ture. With this perspective in sight, it is wise 
in the case of M82X-2 to place constraints on 
the likelihood of the magnetic neutron star pulsar 
systems, so the model of common pulsar will be 
tested critically. 


1.3. Result 


Putting apart the discussion of inherent prob¬ 
lems of the mass scaling of the black holes in 
ULXs, which is beyond the scope of this report, 
we approach the M82X-2 issue from the stand¬ 
point of black holes rather than magnetic neutron 
star pulsar systems. To reconcile the observed un¬ 
usual high pulsed luminosity with the above men¬ 
tioned violation of the Eddington limit, we exam¬ 
ine the physics which is at work in ultraluminous 
pulsations by assuming that M82X-2 is being a 
spinning intermediate-mass black hole (SIMBH). 
If a black hole of intermediate mass will be an ex¬ 
act law of Nature, it i s certainly an attrac tive sce¬ 


nario. Fore example, iKong et al. (1200711 suggest 
that the M82X-2 is a binary system with a black 
hole accretor. Assuming the persistent emission 
is isotropic, the X-ray luminosity ~ 10^°ergs“^ 
implies that the compact object is a > IOOMq 
IMBH in the low/hard state. However, one may 
be tempted to argue truly that most of the issues 
and objections raised above cannot be solved in 
the framework of conventional black hole models, 
to which we refer as phenomenological black hole 
models (PBHMs) (see subsect.2.1). The coherent 
periodicity obviously rules out PBHM, because 1) 
black holes do not radiate; 2) the spinning black 
holes are axisymmetric and have no internal struc¬ 
ture on which to attach a periodic emitter. Orbital 
motion, whether modulating some emission mech¬ 
anism directly or exciting short-period pulsations, 
would decay very quickly due to gravitational radi¬ 
ation. With this in mind, we revisited the MTBH 
which completes PBHM by exploring the most im¬ 
portant processes of spontaneous breaking of grav¬ 
itation gauge symmetry and rearrangement of vac¬ 
uum state at huge energies. One of the purposes 
of this report is to motivate and justify the fur¬ 
ther implications of MTBH framework to circum¬ 
vent the alluded obstacles without the need for 
significant breaking of Eddington limit. We will 
proceed according to the following structure. To 
start with, some latest developments on IMBHs 
are discussed in subsection 1.4. We provide an 
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analysis aimed at clarifying the current situation. 
To make the rest of paper understandable, sec¬ 
tion 2 deals with a brief review of key objectives 
of the MTBH framework, in particular, in rele¬ 
vance to IMBH seeds as a guiding principle. In 
section 3 we set out to examine the ultraluminous 
pulsations powered by M82X-2, constituting mass- 
exchange binary with the 0/B-type donor star. 
We discuss a basic geometry which describes the 
rotating axisymmetric black holes and bring some 
observations which comprise the whole of the case. 
We derive a general profile of pulsed luminosity of 
M82X-2, give a quantitative account of a poten¬ 
tial dynamical mass scaling of M82X-2 and other 
energetics, estimate the mass of companion and 
the orbit parameters of the mass-exchange binary, 
discuss the measured spin-up rate, and calculate 
the torque added to M82X-2 per unit mass of ac¬ 
creted matter. Concluding remarks are presented 
in section 4. 


1.4. Some latest developments on IMBHs 


Transient behavior of M82X-2 with a massive 
companion donor star likely requires an IMBH, 
even though the physical interpretation of the lat¬ 
ter is still controversial. If IMBHs exist they have 
important implications for the dynamics of stellar 
clusters and the formation of supermassive black 
holes. These objects are of particular interest be¬ 
cause this population would fall within the gap be¬ 
tween the two principal black hole populations, i.e. 
super-massive black holes (SMBHs), with masses 
of 10® — 10 ^° Mq, and stellar mass black holes, 
with mas ses of ^ IQMc^. A ROSAT/H RI catalog 
of ULXs ( Colbert fc Mushotzkv 19991) provides a 
clue to fill in the missing population of the IMBHs. 
Currently there are very few w ell studied IMBHs 


two examples are NGC 4395 (IFelioDenko fc Ho 
20031) and POX 52 (iKimth et Ill 9871:iBa.rth" et ai 


20041) . The IMBHs have non-st andard spectra 


which show a cutoff at a few keV (|Stobbart et al. 
20061) . Note that other black holes, stellar as 


well as super-massive, at a few per cent of Ed- 
ding ton luminosity have ha rd power-law-like spec¬ 
tra (IZdziarski et al. 111997^ . Some ULXs show a 
harder, 14 keV thermal compone nt with the corre¬ 
sponding radius of only 304 0 km (|Makishima et al. 
2 OOOI : Stobbart fit al. ll2006h . It is not clear to date 
how IMBHs would form. They have remained ob- 
servationally elusive, with dynamical evidence for 


such objects in large globular clusters still the 
subject of some dispute. On the one hand, to pro¬ 
duce an massive IMBH of 10^ — 10"^ Mq, the core 
collapse of an isolated star in the current epoch is 
not a viable process, which is how the stellar black 
holes are thought to form, because of the lack of 
very massive stars. Their environments lack the 
extreme conditions, i.e. high density and veloc¬ 
ities observed at the centers of galaxies, which 
seemingly lead to the formation of suppermassive 
black holes. There are two conventional scenar¬ 
ios for the formation of IMBHs. The first is the 
merging of stellar mass black holes and other com¬ 
pact objects by means of gravitational radiation. 
The second one is the runaway collision of massive 
stars in dense stellar clusters and the collapse of 
the collision product into an IMBH. But, in fact, 
most ULX host galaxies do not even have stellar 
clusters sufficiently massive and compact to sat¬ 
isfy the requirements for runaway core collapse. 
These objects could be formed inside clusters that 
have since dispersed. However, the evaporation 
timescale of such clusters would be too long to 
explain the observed association of many ULXs 
with young ( < 20 Myr) stellar populations. For 
example, one of the strongest candidates for IMBH 
to date is the hyperluminous X-ray source HLXl 
or 2XMM J011028.1-460421, possibly located in 
the galaxy ES0243-49, or, at least, projected in¬ 
side the = 25.0 mag arcs ec~^ surface bright¬ 
ness isophote of th at galaxy dFarrell et al. I 2009t 
ICodet et al. 20091) . Soria et al. (2010) report a 
discovery of the likely optical counterpart to this 
source. In their analysis of the UV emission in 
ES0243- 49, they determine the optical flux, and 
the X-ray/optical flux ratio, and conclude that 
the X-ray source belongs to ES0243-49, seem¬ 
ingly located inside a massive star cluster. A 
metal-free population HI stars formed in the very 
early Universe could reach masses of a few hun¬ 
dred Mq, above the pair-instabilit y limit, and thus 


may have collapsed into IMBHs (jMadau fc Rees 


200 iD . It was then expected that in young 


and dense star clusters, dynamical friction could 
lead to massive stars sinking towards the center 
and undergoing runaway collisions and mergers 
on timescales 10® yr. Dynamical evidence for 
IMBHs with masses ^ 10"'^ Mq has been pro¬ 
posed for a few globular clusters , e.g., G1 in M31 


at the 1.5cr significance level (IGebhardt et al. 
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20051) . although the issue is still controversial, 


and these results require further extensive and 
careful investigation of the stellar density and 
velocity distributions in globulars, being of the 
highest importance in understanding of the for¬ 


mation of IMBHs (|Anderson. &: van der Marel 


l2010l ). In massive and/or compact globular clus¬ 
ters, a central seed black hole may grow by up 
to a factor of 100 via accretion of gas lost by 
the hrst gen eration of cluster stars in their red- 
giant phase ( Vesperini et al. 2010l) : or, IMBHs 
may wander in the halo of major galaxies, after 
tidal stripping of mergin g satellite dwarfs that 
contained nucley BH s ( King fc Dehnen 1 1200, 51: 
Bellovarv et al. I 2010 l) . Although this seems one 


of the most plausible explanations for the above 
mentioned source HLX-1, in general, currently we 
do not have a complete theoretical interpretation 
of physical nature of growth of black hole seeds. 
Therefore, the MTBH framework may become of 
eminent physical significance in tackling this prob¬ 
lem. 


2. MTBH, Revisited: the implications for 
IMBHs 

The aim of the present section is to recount 
some of the highlights behind the MTBH for the 
benefit of the reader, as a guiding principle to 
make the rest of paper understandable. We shall 
see how a IMBH seed is thought to form. There 
are several important topics not touched upon 
here, which will eventually benefit from a proposed 
gravitation theory. Although some key theoretical 
ideas were introduced with a satisfactory substan¬ 
tiation, we have also attempted to maintain a bal¬ 
ance being overly detailed and overly schematic. 

2.1. Phenomenological model of black 
holes 

From its historical development, up to current 
interests, the efforts in the active galactic nuclei 
(AGN) physics have evoked the study of a ma¬ 
jor unsolved problem of how efficiently such huge 
energies observed can be generated. This energy 
scale severely challenges conventional source mod¬ 
els. The huge energy release from compact regions 
of AGNs requires extremely high efficiency (typi¬ 
cally > 10 per cent) of conversion of rest mass to 
other forms of energy. This led early on to sug¬ 


gestions that AGNs are powered by SMBHs, with 
masses of 10® —10^° Mq . A complex study of AGN 
evolution requires a comprehensive understanding 
of a new field of astrophysics, so-called black hole 
demography. This studies the important phe¬ 
nomena of the birth and growth of black holes, 
the merging of stellar mass black holes and their 
evolutionary connection to other objects in the 
Universe. These ideas gather support especially 
from a breakthrough made in recent observational 
and computational efforts on understanding of 
coevolution of black holes and their host galax¬ 
ies, particularly through self-regulated growth 
and feedback from accret i on-pow e red outflows 


(e.g. 

Nataraian &: Treister ( 

2009): Vestergaard 

(2004 

): Volonteri et oi. (2008 

):IVolonteri & Nataraian 

(2009 

): Volonteri (2010l): IShankar et al. 1 (2009): 

Kellv et al. (20101): Na 

taraian (2011): Treister & Metran 

(2012 

): Willott et at. ( 

2010r):lDevecchi & Volonteri 1 

(2009 

). In these models, at early times the prop- 


erties of the assembling SMBH seeds are more 
tightly coupled to properties of the dark matter 
halo as their growth is driven by the merger his¬ 
tory of halos. 

General relativity (GR) has stood the test 
of time and can claim remarkable success, al¬ 
though there are serious problems concerning the 
role of singularities or black holes. This state 
of affairs has not much changed up to present 
and proposed abundant models are not conduc¬ 
tive to provide non-artificial and unique recipe 
for resolving controversial problems of energy- 
momentum conservation laws of gravitational in¬ 
teracting fields, localization of energy of gravita¬ 
tion waves and also severe problems involved in 
quantum gravity. A tacit assumption of theoreti- 


cal interpretation of aforementioned astrop 
scenarios is a general belief (e.g. Rees ( 

ivsic 

1984 

al 

); 

); 

Shaniro (2005) 

Munvaneza & Biermann 1 ( 

2006 

Alexander & Nataraian (|2014|) and references 


therein) reinforced by statements in textbooks, 
that a longstanding standard PBHM can describe 
the growth of accreting black hole seeds. In the 
framework of GR, the PBHM implies the most 
general Kerr-Newman black hole model, with the 
integral parameters of total mass (M), angular 
momentum (J) and charge (Q), still has to put 
in by hand. But such beliefs are suspect and 
should be critically re-examined. Even though 
being among the most significant advances in as- 
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trophysics, it is rather surprising that PBHM is 
routinely used to explore this problem as though 
it cannot be accepted as convincing model for ad¬ 
dressing the black hole growth, because in this 
framework the very source of gravitational field 
of the black hole is a kind of meaningless cur¬ 
vature singularity at the center of the stationary 
black hole, which is hidden behind the event hori¬ 
zon. The theory breaks down inside the event 
horizon which is causally disconnected from the 
exterior world. Either the Kruskal continuation of 
the Schwarzschild (J = 0, Q = 0) metric, or the 
Kerr {Q = 0) metric, or the Reissner-Nordstrom 
(J = 0) metric, show that the static observers 
fail to exist inside the horizon. Any object that 
collapses to form a black hole will go on to col¬ 
lapse to a singularity inside the black hole. Any 
timelike worldline must strike the central singu¬ 
larity which wholly absorbs the infalling matter. 
Therefore, the ultimate fate of collapsing matter 
once it has crossed the black hole surface is un¬ 
known. So, one should deliberately forebear from 
presumption of exotic hypothetical growth behav¬ 
ior of black holes, which seems nowhere near true 
if one applies the phenomenological model. This 
ultimately disables any accumulation of matter in 
the central part and, thus, neither the growth of 
black holes nor the increase of their mass-energy 
density could occur at accretion of outside mat¬ 
ter, or by means of merger processes. Yet it is 
still thought provoking how can one be sure that 
some hitherto unknown source of pressure does 
not become important at huge energies and halt 
the infinite collapse. That in PBHM there is no 
provision for growth behavior of black holes, is 
because in it one assigned only an insufficient at¬ 
tributes to this. The PBHM is a rather restricted 
model and one needs to realise that if one can 
gain insight into exploring a new process of spon¬ 
taneous breaking of gravitation at huge energies, 
one has then made room for growth and merging 
behavior of black holes. To fill the void which the 
standard PBHM presents, thus, one plausible idea 
to innovate the solution to alluded key problems 
would appear to be such breaking mechanism, and 
thereof for that of rearrangement of vacuum state 
itself. 


2.2. Proposed gravitation theory 


In this subsection we present a brief outline 
of the key ideas of the underlying theory, which 
yields the gravitational interaction at huge en¬ 
ergies drastically different from earlier suggested 
schemes. This theory involves a drastic revision 
of a role of local internal symmetries in physi¬ 
cal concept of curved geometry, and explores the 
most important processes of spontaneous break¬ 
ing of gravitation gauge symmetry and rearrange¬ 
ment of v acuum state. It was originally pro¬ 
posed by ( Ter-Kazarian I ( 19971 . 12001 . 2010l l and 
references therein), and th oroughly discussed in 
I Ter-Kazarian 2014 . 2015ah . This extension, suit¬ 
able for applications in extremely high energy as¬ 
trophysics, is a bold assumption in its own right. 
In its present formulation, this theory exploits 
the language of the fundamental geometric struc¬ 
ture- distortion gauge induced fiber-hundle, which 
provides a modified gravitational theory, as a 
corollary of the spacetime deformation/d is tortion 
framewo rk, suggested b y Ter-K azarian I (2015a) 
(see also Ter-Kazarian I (1201 iL 1201 bbl lTlt should 
be emphasized, that the standard Riemannian 
(and its extensions) space interacting quantum 
field theory cannot be a satisfactory ground for 
addressing the problems in quest. The difficul¬ 
ties associated with this step are notorious, how¬ 
ever, these difficulties are technical. In the main, 
they stem from the fact that Riemannian ge¬ 
ometry, in general, does not admit a group of 
isometries, and that it is impossible to define 
energy-momentum as Noether local currents re¬ 
lated to exact symmetries. This, in turn, posed 
severe problem of non-uniqueness of the physi¬ 
cal vacuum and the associated Fock space. A 
definition of positive frequency modes cannot, in 
general, be unambiguously fixed in the past and 
future, which leads to \in \out >, because 
the state \in > is unstable against decay into 
many particle \out > states due to interaction 
processes allowed by lack of Poincare invariance. 
A non-trivial Bogolubov transformation between 
past and future positive frequency modes implies 
that particles are created from the vacuum and 
this is one of the reasons for \in \out >. 
Keeping in mind aforesaid, we develop on the 
framework of the general gauge principle, as 
from first principles'. We consider the principal 
fiber bundle with the semi-Riemannian spacetime. 
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14 , as the base space, and with the structure 
group, Gy, generated by the hidden local internal 
symmetry of two-parameter abelian local group 
t/'°=( 2 ) = U{\)y X 17(1) = U{1)y X diag[SU{2)]. 
This group is implemented on the flat space, M 4 . 
The collection of matter fields of arbitrary spins 
take values in standard fiber. To involve a drastic 
revision of the role of gauge fields in the physical 
concept of the spacetime deformation/distortion, 
we generalize the standard gauge scheme via the 
concept of distortion gauge field which acts on 
the external spacetime groups. The group 17*°°(2) 
has two generators, the third component of 
isospin T related to the Pauli spin matrix and 
hypercharge Y implying Q'* = where 

is the distortion charge operator assigning the 
number -1 to particles, but -1-1 to anti-particles. 
We connect the structure group Gy, further, to 
the nonlinear realization of the Lie group G n 


of dis t ortion of the spacetime le.g iTer-Kazarian 
(12014 [ioili)), i.e. we extend the curvature of 
the spacetime continuum to general distortion as 
the theory of spontaneous breaking of distortion 
symmetry. This can be achieved by non-linear 
realizations of the distortion group Gd ■ The non¬ 


nomenological Lagrangians (IColeman et al. 

1967; 

Gallan et al. 

Il969l: Salam & Strathdee 

196S; 

Isham et al. 

I 971 I: Ogievetskv 20091: Volkov 

2010; Weinberg 2010ll provides a wav to deter- 


mine the transformation properties of fields de¬ 
fined on the quotient space. In accord, we treat 
the distortion group Gd and its stationary sub¬ 
group H = SO{3), respectively, as the dynamical 
group and its algebraic subgroup. The fundamen¬ 
tal field is distortion gauge field and, thus, all 
the fundamental gravitational structures in fact - 
the metric as much as the coframes and connec¬ 
tions - acquire a distortion-gauge induced theo¬ 
retical interpretation. We study the geometrical 
structure of the space of parameters in terms of 
Cartan’s calculus of exterior forms and derive the 
Maurer-Cartan structure equations. There is a 
reciprocity, treated in the Maurer-Cartan nonlin¬ 
ear structure equations, between the formal math¬ 
ematical structure of distortion fields on the one 
hand, and the Goldstone fields, on the other. The 
group {7*°°(2) entails two neutral gauge bosons of 
C/(l), or that coupled to T^, and of 17(l)y, or that 
coupled to the hypercharge Y. Spontaneous sym¬ 


metry breaking can be achieved by introducing 
the neutral complex scalar Higgs field. Minimiza¬ 
tion of the vacuum energy fixes the non-vanishing 
vacuum expectation value, which spontaneously 
breaks the theory, leaving the U{l)d subgroup 
intact, i.e. leaving one Goldstone boson. Conse¬ 
quently, the left Goldstone boson is gauged away 
from the scalar sector, but it essentially reappears 
in the gauge sector providing the longitudinally 
polarized spin state of one of gauge bosons that 
acquires mass through its coupling to Higgs scalar. 
Thus, the two neutral gauge bosons were mixed to 
form two physical orthogonal states of the mass¬ 
less component of distortion field, (a) {Ma = 0 ), 
which is responsible for gravitational interactions, 
and its massive component, (a) (Mg ^ 0 ), which 
is responsible for the, so-called inner distortion 
(ID)-regime. Hence, a substantial change of the 
properties of the spacetime continuum besides the 
curvature may arise at huge energies. This theory 
is renormalizable, because gauge invariance gives 
conservation of charge, also ensures the cancela¬ 
tion of quantum corrections that would otherwise 
result in infinitely large amplitudes. Without care¬ 
ful thought we expect that in this framework the 
renormalizability of the theory will not be spoiled 
in curved space-time as well, because, the infini¬ 
ties arise from ultra-violet properties of Feynman 
integrals in momentum space which, in coordinate 
space, are short distance properties, and locally 
(over short distances) all the curved spaces look 
like maximally symmetric (flat) space. 

For brevity, we will not be further concerned 
with the actual details of this comprehensive theo¬ 
retical framework, but only use it as a backdrop to 
validate the theory with more observational tests. 
For details, the interested reader is invited to con¬ 
sult the original papers. Discussed gravitation the¬ 
ory is consistent with GR up to the limit of neu¬ 
tron stars. But MTBH manifests its virtues ap¬ 
plied to the physics at huge energies. Whereas a 
significant change of properties of spacetime con¬ 
tinuum, ID-regime (at a ^ 0), arises simultane¬ 
ously with the strong gravity (a 7 ^ 0). The Gomp- 
ton length of the ID-field d{r) is Ag = — 

0.4/m. 






































2.3. Crossing event horizon from inside 
the black hole 


(ao) a nd ID (a) fields can b e obtained in Feynman 
gauge (jTer-Kazarian Il200lh as 


Let us to discuss in more details some principle 
issues in use. Recall that according to GR, the 
event horizon is impenetrable barrier for crossing 
from inside the black hole because of a singular¬ 
ity arisen at the Schwarzschild radius (radius of a 
non-rotating black hole) Rg = 2GM/c^ = 2.95 x 
cm. But in the framework of MTBH, 
this barrier disappears when a matter found in 
the ID-region of space-time continuum has under¬ 
gone phase transition of II-kind and, thus, it trans¬ 
formed to proto-matter. To obtain some feeling 
about this phenomena, below we give more de¬ 
tailed explanation. In the framework of MTBH, 
we are led to construct a formalism of unitary 
mapping of the fields and their dynamics from M 4 
to V 4 , and vice versa. The field equations follow 
at once from the total gauge invariant Lagrangian 
in terms of Euler-Lagrange va riations, respectively 
on bo th 14 and M 4 spaces ( Ter-Kazarian I 20l4 


201 bal l. The Lagrangian of distortion gauge field 


defined on the flat space is undegenerated Killing 
form on the Lie algebra of the group in 

adjoint representation, which yields the equation 
of distortion field. In the reminder of this sub¬ 
section the quantities referred to 14 are denoted 
by wiggles, and left without wiggles if they cor¬ 
respond to M 4 . We are interested in the case of 
a spherical-symmetric gravitational field oo(r) in 
presence of one-dimensional space-like ID-held a. 
In the case at hand, one has the group of mo¬ 
tions SO{3) with 2D space-like orbits where 
the standard coordinates are 6 and The sta¬ 
tionary subgroup of SO (3) acts isotropically upon 
the tangent space at the point of £phere 5^ of ra¬ 
dius r. So, the bundle p ■. V 4 ^ B? has the hber 
= p~^{x), X G V 4 with a trivial connection 
on it, where R^ is the quotient-space V 4 /SO{ 3 ). 
Considering the equilibrium conhgurations of de¬ 
generate barionic matter, we assume an absence of 
transversal stresses and the transference of masses 
in 14 


Tl = Ti = Ti = -P(?), T° = -p(?), ( 1 ) 

where P{r) and p(r) (r G R^) are taken to de¬ 
note the internal pressure and macroscopic den¬ 
sity of energy dehned in proper frame of reference 
that is being used. The equations of gravitation 


Aao = ^{aoo^pir}- 


933 Qao 


■ Ti 

933 Qao 


■ Ti ^ 

1^22 Qag 


P 


(^} 


(A-A„ 2 )a= i{goo^p(? 0 - 


( 2 ) 


933-§^ + gii-§r + 922-§r 


P(?)} 


where a diffeomorphism r(r) : M 4 —>■ 14 is given 
r = r — Rg/A. A distortion of the basis e in 
the ID regime, in turn, yields the transforma¬ 
tions of Poincare generators of translations. Given 
an explicit form of distorted basis vectors, it is 
straightforward to derive the laws of phase transi¬ 
tion for individual particle found in the ID-region 
(xo = seao = 0 , x = sea ^ 0 ) of the space-time 
continuum, where a coupling constant ae relates to 
Newton gravitational constant G as se^ = SttG/c^. 
The Poincare generators P^ of tr anslations un¬ 
dergone follo wing transformations ( Ter-Kazarian I 
20niU2015ah : 


E = E, P 1.2 = Pi .2 cos 63 


tan 03 me, m = 


(m - tan 63 


P 3 = P 3 - 
2 


(3) 


sin^ 63 ■ 


— tan^ 03 




where E, P, m and E, P, fh are ordinary and 
distorted energy, momentum and mass at rest, 
tan 03 = —X. Hence the matter found in the ID- 
region (a ^ 0) of space-time continuum has under¬ 
gone phase transition of H-kind, i.e., each particle 
goes off from the mass shell- a shift of mass and 
energy-momentum spectra occurs upwards along 
the energy scale. The matter in this state is called 
proto-matter with the thermodynamics differed 
drastically from the thermodynamics of ordinary 
compressed matter. The resulting deformed met¬ 
ric on 14 in holonomic coordinate basis takes the 
form 

000 = ( 1 -a^o)^+^^ 9g.v=^ 


033 = - [(1+ a;o)^ , 011 =-r^, (4) 

922 = -r^ sin^ 0. 
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The metric o clearly shows that: a singularity 
at intersection of proto-matter disk with the event 
horizon (at xq = 1 ) disappears where a massive 
component of ID-field is not zero (x ^ 0), and 
hence the crossing event horizon from inside of 
black hole at such conditions is allowed. 


2.4. Rational 


Consisting of the proto-matter core and the 
outer layers of ordinary matter, the SPC config¬ 
uration is the spherical-symmetric distribution of 
barionic matter in many-phase stratified states. 
A layering is a consequence of the onset of dif¬ 
ferent regimes in equation of state. Each con¬ 
figuration is defined by the two free parameters 
of central values of particle concentration n(0) 
and dimensionless potential of space-like ID-held 
a;(0). The interior gravitational potential x(f*^{r) 
matches into the exterior one x(f*(r) at the sur¬ 
face of the conhguration. The central value of 
the gravitational potential a:o(0) can be found by 
reiterating integrations when the sewing condi¬ 
tion of the interior and exterior potentials holds. 
The simulations conhrm in brief the following sce¬ 
nario. The energy density and internal pressure 
have sharply increased in proto-matter core, with 
respect to corresponding central values of neutron 
star, proportional to gravitational forces of com¬ 
pression. This counteracts the collapse and equi¬ 
librium holds even for the masses up to ~ 10 ^°Mq. 
Encapsulated in a complete set of equations of 
SPC-conhguration, the SPC is a robust structure 
that has stood the tests of the most rigorous the¬ 
oretical scru tinies of its stabilitv, which was a cen¬ 
tral issue in (|Ter-Kazarian et al. II2007I) . Minimiz¬ 
ing the total energy gives the equilibrium config¬ 
urations. The second derivative of total energy 
gives stability information. Although a relativ¬ 
ity tends to destabilize configurations, however, a 
numerical integrations of the stability equations 
of SPC clearly proves the stability of resulting 
cores. Due to it, the stable equilibrium holds as 
well in outward layers and, thus, an accumula¬ 
tion of matter is now allowed around the stable 
SPC. To emphasize the distinction between phe¬ 
nomenological and microscopic black hole models, 
we present their schematic plots in Fig.[Tl to guide 
the eye. One of the most remarkable drawback of 
MTBH is that the central singularity cannot oc¬ 
cur, which is replaced by finite though unbeliev- 


PBHM MTBM 




Fig. 1.— Left panel. Phenomenological model of 
non-spinning black hole. The meaningless singu¬ 
larity occurs at the center inside the black hole. 
Right panel: Microscopic model of non-spinning 
black hole, with the central stable SPC. An in¬ 
falling matter with the time forms PD around the 
SPC. In final stage of growth, a PD has reached 
out the edge of the event horizon. Whereas a 
metric singularity inevitably disappears and UHE 
neutrinos may escape from event horizon to out¬ 
side world through vista, i.e. a thin belt area 
S = 2-KRgd of intersection of proto-matter disk 
with the event horizon, with opening angle Oi,. 
Accepted notations: EH=Event Horizon, 

AD=Accretion Disk, SPC=Superdense Proto¬ 
matter Core, PD=Proto-matter Disk. 


ably extreme conditions held in the stable SPC 
where, nevertheless, static observers are existed. 
The SPC is always found inside the event hori¬ 
zon sphere, therefore it could be observed only 
in presence of accreting matter. The SPC sur¬ 
rounded by the accretion disk presents the micro¬ 
scopic model of AGN. The SPC accommodates 
the highest energy scale up to hundreds ZeV in 
central proto-matter core of suppermasive black 
holes, which accounts for the spectral distribu¬ 
tion of the resulting radiation of galactic nuclei. 
An external physics of accretion onto the black 
hole in the earlier half of its lifetime is identical to 
the processes in Schwarzschild’s model. A crucial 
difference in the model context between the phe¬ 
nomenological and microscopic black holes comes 
in when one looks for the spontaneous breaking 
of gravitation gauge symmetry at huge energies, 
and thereof making room for growth and merg¬ 
ing behavior of black holes. We have shown that 
the AGN evolution and black hole growth, when 
properly analysed, agree with our claims. 
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2.5. The intermediate mass, pre-radiation 
time and initial redshift of seed black 
hole 


We argue that in the framework of MTBH, the 
black hole seeds might grow driven by the ac¬ 
cretion of outside matter when they were getting 
most of their masses. An infalling matter with 
time forms a thin proto-matter disk around the 
proto-matter core tapering off faster at reaching 
out the event horizon. As a a metric singularity 
inevitably disappears and the ZeV-neutrinos, pro¬ 
duced via simple or modified URCA0 processes in 
deep layers of SPC and proto-matter disk, may es¬ 
cape from event horizon to outside world through 
a thin belt area S = 27r Rgd, even after the strong 
neutrino trapping. The thickness of proto-matter 
disk at the edge of event horizon is denoted by d. 
The neutrinos are collimated in very small open¬ 
ing angle 9^ — Sd = ^ 1. The trapping is 

due to the fact that as the neutrinos are formed in 
proto-matter medium at super-high densities they 
experience greater difficulty escaping from it be¬ 
fore being dragged along with the matter, i.e. the 
neutrinos are trapped comove with matter. In this 
framework we introduced a notion of pre-radiation 
time (PRT) of black hole which is referred to as a 
lapse of time Tbh from the birth of black hole till 
neutrino radiation - provision for the earlier half 
of the lifetime of black hole: 


rr. Md 

^BH = -r-- 

M 


(5) 


Here Md is the mass of proto-matter disk, M is 
an appropriately averaged mass accretion rate. To 
render our discussion here a bit more transparent, 
we present the relation between typical PRT ver¬ 
sus bolometric luminosity of suppermassive black 
holes as follows: 


Tbh ^ 0.32— 
rov 


Mbh V 

Mq J Lhoi ^ 


( 6 ) 


where Rd is the radius of the proto-matter core, 
fov = 13.68km. At times > Tbh, the black 
hole no longer holds as a region of spacetime 


^G. Gamow was inspired to name the process URCA after 
the name of a casino in Rio de Janeiro, when M. Schenberg 
remarked to him that ’’the energy disappears in the nucleus 
of the supernova as quickly as the money disappeared at 
that roulette table” 


that cannot communicate with the external Uni¬ 
verse. In this framework, we computed the fluxes 
of ZeV-neutrinos from plausible accreting super- 
massive Un£k_holesclosely l inking with the 377 
AGNs ( Ter-Kazarian I 2014 . 2015all . In accord, 
the AGNs are favored as promising pure neutrino 
sources because the computed neutrino fluxes are 
highly beamed along the plane of accretion disk, 
and peaked at high energies and collimated in 
very small opening angle 0;, ~ "C 1- While 
hard to detect, the extragalactic ZeV-neutrinos 
may reveal clues on the puzzle of the origin of 
ultra-high energy cosmic rays, as they have the 
advantage of representing unique fingerprints of 
hadron interactions and, therefore, can initiate 
the cascades of UHE-particles with energies ex¬ 
ceeding 1.0 X 10^ ° eV (comprehensive r eview s 


can be found in |CMtellina_^E)qnatcL 

Letessier-Selvon fc Staney (12011 ): Sid 


Kotera fc Olinto I (l20inll : ISemikoz I (l20inlP . 


( 2012 ); 

dmH); 


Consequently, we have studied a growth of 
proto-mat ter disk and derived the mass of black 
hole seed ( Ter-Kazarian |[^15a ) - provision for the 
second half of the lifetime of black hole: 


Mr. 


M, 


BH 


Mr. 


Rd Mbh \ 

1 _ 2.305 ) , (7) 

rov Mq 


and initial redshift 


^Seed 


z + Ho Tbh, 


( 8 ) 


where Hq is the Hubble’s constant. Whereas inter¬ 
preting the redshift as a cosmological Doppler ef¬ 
fect, and that the Hubble law could most easily be 
understood in terms of expansion of the Universe, 
as a further consistency check, we then turned the 
problem around by asking the purely academic 
question of principle what could be the initial red¬ 
shift, of seed black hole if the mass, the 

luminosity and the redshift, z, of black hole at 
present time are known. Equation (0 holds at 
< 0.023-^, where i?® is the gravitational 
radius of black hole seed. We have undertaken a 
large series of numerical simulations with the goal 
to trace an evolution of the mass assembly history 
of 377 accreting supermassive black hole seeds in 
AGNs to the present time and examine the observ¬ 
able signatures today. Given the redshifts, masses 
and luminosities of these black holes at present 
time collected from the literature, we compute the 
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initial redshifts and masses of the corresponding 
black hole seeds. 

Having gained some insight into the supermas- 
sive black hole physics, let us now comment briefly 
on the implications for IMBHs. Equation © 
shows that the seed of IMBH can be formed in 
the stellar mass black hole region if the condition 


Rd 

rov 


< 0.434 


Mq 

Mbh 



(9) 


holds. A creation of IMBHs may be possible in the 
young stellar populations that we generally find 
ULXs co-habiting with. The IMBH with mass, 
say Mbh — 1-0 x 10^ Mq, may originate from the 
seed with stellar mass of ^ I Mq at Rd/rov — 
4.3 X 10“"^, and will reach to the finale state after 
a lapse of growing time of order Tbh — Mbh/M. 


3. The mass-exchange binary containing 
M82X-2 


As in the case of neutron stars, we expect 
that accreting black holes are fast spinning ob¬ 
jects. For the self-contained arguments, we need 
to extend the preceding algorithm of non-spinning 
MBHM to its spinning counterpart, which is 
almost a matter of routine, to change the ge¬ 
ometry of static SPC to a more general one, 
describing axisymmetric rotating SPC. For the 
standard calculations of rapidly rotating rela- 


fer to 

fe.er. Carter ( 

1970); 

Bardeen ( 

1970 

); 

); 

Komatsu et al. I(ll989ll: 

Cook et al I(ll992. 

1994 

Bonazzola et al. 1 (ll993h;lBonazzola & Gorsoulhon 

(I 994 II: 

Sterffioulas | ( 2003f)l and references therein. 


3.1. The geometry of rotating 8Lxisymmet- 
ric SPC 

In this subsection we will collect together the 
results which are required later. The non-spinning 
SPC is static and spherically symmetric. So, 
one needs to be clear about more general ge¬ 
ometry which can describe rotating axisymmet¬ 
ric SPCs. The principle foundation of the spin¬ 
ning configurations first comprises the following 
additional distinctive features with respect to non¬ 
spinning ones: I) Rapid rotation causes the shape 
of the SPC to be flattened by centrifugal forces- 
flattened at poles and bulged at equator (oblate 


spheroid, which is second order effect in the rota¬ 
tion rate). 2) A rotating massive SPC drags space 
and time around with it. The local inertial frames 
are dragged by the rotation of the gravitational 
field, i.e. a gyroscope orbiting near the SPC will 
be dragged along with the rapidly rotating SPC. 
This is probably the most remarkable feature that 
could serve as a lin k with the general desc ription of 
spacetime (also see Ter-Kazarian I ( 20I2[B . Beside 
the geodetic procession, a spin of the body pro¬ 
duces in addition the Lense-Thirring procession. 
To look into the future, measurement of the gyro- 
gravitational r atio of par ticle would be a further 
step, see e.g. (iNi 11201011 and references therein, 
towards probing the microscopic origin of gravity. 
Let the world coordinate t{= x^) be the time (in 
units of c), and (/){= x^) be the azimuthal angle 
about the axis of symmetry. Moreover, a metric 
of a two-space with positive or negative signature, 
{x‘^,x^), can always be brought to the diagonal 
form by a more coordinate transformation. Since 
the source of gravitational field has motions that 
are pure rotational about the axis of symmetry, 
then the energy-momentum tensor as the source of 
the metric will have the same symmetry. Namely, 
the space would be invariant against simul¬ 
taneous inversion of time t and azimuthal angle 
(f). The 3 -|- 1 formalism is the most commonly 
used approach in which, as usual, spacetime is de¬ 
composed into the one parameter family of space¬ 
like slices- the hypersurfaces E^. The study of a 
dragging effect is assisted by incorporating with 
the soldering tools in order to relate local Lorentz 
symmetry to curved spacetime. These are the lin¬ 
ear frames and forms in tangent fiber-bundles to 
the external general smooth differential manifold, 
whose components are so-called tetrad (vierbein) 
fields. Whereas, the has at each point a tan¬ 
gent space, TxM-i, spanned by the anholonomic 
orthonormal frame field, e, as a shorthand for 
the collection of the 4-tuplet (eo = exp(—:/) {dt + 
ujd^), ei = exp{-^)d^, ei = exp(-^2)i92, 63 = 
exp(—/rs) 93 ), where ta = is the the 

soldering form between the tangent space and the 
spacetime manifold. This is called a Bardeen ob¬ 
server, locally nonrotating observer, or the local 
Zero Angular-Momentum Observers (ZAMO), i.e. 
observers whose worldlines are normal to the hy¬ 
persurfaces defined by constant coordinate time, 
t = const, also called Eulerian observers. Here we 
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use Greek alphabet (/x, i/, p,... = 0,1,2,3) to de¬ 
note the holonomic world indices related to A^ 4 , 
and the first half of Latin alphabet (a, &, c,... = 
0,1,2,3) to denote the anholonomic indices re¬ 
lated to the tangent space. The frame field, e, 
then defines a dual vector, -d, of differential forms. 


-d = 


\ 


, as a shorthand 


/ = exp 1 / dt 

'd^ = exp tj} {dip — ojdt) 

= exp p 2 dx^ 

/i 3 dx^ J 

for the collection of the d^ = ,, dx^, whose val- 

ues at every point form the dual basis, such that 
ea\d^ = 6 ^, where J denoting the interior prod¬ 
uct, namely, this is a C'°°-bilinear map J : 11^ —>■ 
with denotes the C^-modulo of differen¬ 

tial p-forms on Ai 4 . In components ^ = S^. 
The norm ds of infinitesimal displacement dx^ on 
Ad 4 , describing the stationary and axisymmetric 
spacetimes, reads 


ds : = ed = ef_i ^ d^ G M 4 . ( 10 ) 

Therefore, the holonomic metric on the space Ad 4 
can be recast into the form 


9 = g^.d>^<^d''= g{ef„e,)d^^d-', ( 11 ) 

with the components g^^i, = g{e^, e^) in dual holo¬ 
nomic basis = dx^}. That is 

ds^ = — exp{2v)dt^ + exp{2ip){d(p — ujdt)‘^+ 

exp{2gL2){dx'^Y + exp(2/i3)(da;^)^, 

( 12 ) 

where the five quantities v,ip,u},g ,2 and /X 3 are 
only functions of the coordinate x^ and x^. In 
the case at hand, the metric function w is the an¬ 
gular velocity of the local ZAMO with respect to 
an observer at rest at infinity. Thereby the red- 
shift factor a = exp v is the time dilation factor 
between the proper time of the local ZAMO and 
coordinate time t along a radial coordinate line, 
i.e. the redshift factor for the time-slicing of a 
spacetime. In accord, all the geometrical objects 
are split into corresponding components with re¬ 
spect to this time-slice of spacetime. In partic¬ 
ular, the splitting of manifold Ad 4 into a folia¬ 
tion of three-surfaces will induce a corresponding 
splitting of the affine connection, curvature and, 
thus, of the energy-momentum tensor which can 
as well be written in terms of the energy density 
E = T{n,n) = measured by an adapted 


Eulerian observer of four-velocity n^, the momen¬ 
tum flow Jot = and the correspond¬ 
ing stress tensor Sap = (S = 5'“), that 

is T°‘d = En°‘nl^ + + J°‘nd + 5'“^. Here 

is the unit orthogonal vector to the hypersur¬ 
face Ej, whereas the spacetime metric g induces 
a first fundamental form with the spatial met¬ 
ric jotp on each E* as 'jap = gap + naUp. The 
form (13 includes one gauge freedom for the co¬ 
ordinate choice. For the spherical type coordi¬ 
nates x^ = 9 and x^ = r, for example, so-called 
quasi-isotropic gauge corresponds to = 0 and 
'ygg = r'^jrr- Then, one may define the second 
fundamental form which associates with each vec¬ 
tor tangent to E*, and the extrinsic curvature of 
the hypersurface Et as minus the second funda¬ 
mental form. Aftermath, one can define the usual 
Lorentz factor W = —n^u'^ = avP for a per¬ 
fect fluid with conventional stress-energy tensor 
— {pd~ P)u^u'' + Pg^'', where p is the energy 
density and P is the pressure. Hence E = W^{p + 
P) — P and J* = {E + P)v'‘, where the fluid three- 
velocity u*(i = 1, 2, 3) implies u* = LE(u® — pE j a). 
Thereby the resulting stress tensor can be written 
Sij = {E + P)v^Vj + P^ij. The four-velocity for 
rotating fluid reads u = u^{d/dt) + where 

H is the fluid angular velocity as seen 

by an inertial observer at rest at infinity. It is 
convenient to give the equations in the isotropic 
gauge, exp( 2 /X 2 ) = exp( 2 /i 3 ) = exp(2^), in the 
cylindrical coordinates dg = dx'^, dz = dx^, or 
in pseudospherical coordinates r d9 = dx^ and 
dr = dx^, where the cylindrical radius is writ¬ 
ten exp{ip) = r sm9 B exp(—jx), with B denoting 
a function of r and 6 only. Therefore, the met¬ 
ric m becomes 

ds^ = — exp{2i')dt^ + B^r^ sin^ 9 exp{—2h'){d(p— 

ujdtY + exp{2p){dr‘^ + r‘^d9'^). 

(13) 

Consequently, the components of the energy - 
momentum tensor of matter with total den¬ 
sity p and pressure P are given in the non¬ 
rotating anholonomic orthonormal frame as T°'^ = 
= W^{p + PV^), T^^=W^{p + 

PV2)^ pOl ^ p^y p22 ^ p33 ^ p^ 

with its trace T = —p + 3P, where V is the ve¬ 
locity (in units of c) with respect to the Bardeen 
observer V = ^ so IE = —^==^. However, 

at this point we cut short and, in what follows. 
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Perpendicular 
to Orbit 



Fig. 2.— A schematic SIMBH model of M82X-2 constituting mass-exchange binary with the 0/B-type donor 
star. The angle i is the binary inclination with respect to the plane of the sky. No eclipse condition holds. In 
final stage of growth, PD has reached out the edge of the EH. The thermal defuse blackbody X-rays beams 
may escape from SIMBH through a thin belt area S = 2'KRgd to outside world that sweep past Earth like 
lighthouse beacons. Parameters of a binary system is viewed in the orbital plane. The picture is not to scale. 
Accepted notations: EH=Event Horizon, SPC=Superdense Proto-matter Core, PD=Proto-matter Disk. 


we will refrain from providing further lengthy de¬ 
tails of the mathematical apparatus of proposed 
gravitation theory at huge energies and rigorous 
solution of the extended equations describing the 
spinning MBHM, which is beyond the scope of 
this report. This entire set includes the gravi¬ 
tational and ID field equations of the time-slices 
generated by the above stress-energy tensor, the 
angular momentum equation with the momen¬ 
tum flow determining the frame-dragging poten¬ 
tial, the hydrostatic equilibrium equation and the 
state equation specified for each domain of many 
layered spinning SPC- configurations. Theoretical 
evolutionary paths of this type which are suit¬ 
able for rigorous comparison with the behavior of 
observed sources will be separate topic of com¬ 
prehensive investigation elsewhere. But some ev¬ 
idence for a simplified physical picture, without 
loss of generality, is highlighted in the rest of this 
section where we extend preceding developments 
of MTBH in concise form, without going into the 
subtleties, as applied to the initially rigid-body 
spinning IMBH configuration of angular velocity 
D. 


3.2. The X-ray radiation from SIMBH 

If this were the case, eventually the spinning 
proto-matter core and a thin co-spinning proto¬ 
matter disk driven by accretion would be formed. 
The evolution and structure of a proto-matter disk 
is largely determined by internal friction. Be¬ 
fore tempting to build a physical model in quest, 
the other features of SIMBH configuration also 
need to be accounted which comprise the whole 
of the case. The fact that the rotational energy 
has a steeper dependence on the radius of the 
compact object than the internal energy in the 
relativistic limit is quite significant. Equilibrium 
can always be achieved for massive configurations 
with nonzero angular momentum by decreasing 
its radius. Also, there are two characteristic fea¬ 
tures that distinguish a spinning relativistic SPC- 
conhguration from its non-spinning counterpart: 
3) The geodetic effect, as in case of a gyroscope, 
leads an accretion stream to a tilting of its spin 
axis in the plain of the orbit. Hence a proto-matter 
disk will be tilted from the plane of accretion on 
a definite angle 5 towards the equator. 4) Besides 
the UHE neutrinos, produced in the deep interior 
layers of superdense proto-matter medium as in 
case of non-spinning model, the additional ther- 
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mal defuse blackbody radiation is released from 
the outer surface layers of ordinary matter of spin¬ 
ning SPC and co-spinning proto-matter thin disk. 
All of the rotational kinetic energy is dissipated 
as thermal blackbody radiation. This is due to 
the physical condition that these layers optically 
thick and, eventually, in earlier half of the life¬ 
time of spinning black hole, at times < Tbhi the 
strict thermodynamic equilibrium prevails for this 
radiation because there would be no net flux to 
outside of event horizon in any direction. That is, 
the emission from the isothermal, optically thick 
outer layers at surface is blackbody, which is the 
most efficient radiation mechanism. This radia¬ 
tion is free of trapping. With this guidelines to 
follow, we may proceed in relatively simple way 
toward first look at some of the associated physics 
and can be quick to estimate the physical charac¬ 
teristics of mass-exchange X-ray binaries. Exam¬ 
ining the pulsations revealed from M82X-2, as a 
working model we assume the source of the flashes 
to be a self-gravitating SIMBH resided in the fi¬ 
nal stage of growth. This implies that a thin co¬ 
spinning proto-matter disk has reached out the 
edge of the event horizon. Due to it a metric 
singularity inevitably disappears at the intersec¬ 
tion of proto-matter disk with the event horizon. 
Hence, the crossing event horizon at such condi¬ 
tions from inside of black hole is allowed. This 
general behavior is very robust and depends very 
little on the details of the model of SPC. A prin¬ 
ciple physical properties of this phenomena for 
non-rotating SPC are already discussed in subsec¬ 
tion 2.3. Without being carefully treated, even 
though these properties are of great significant for 
a rotating SPCs too. Actually, by virtue of sub¬ 
sect. 3.1, they will retain for a more general ge¬ 
ometry if metric ([4]) will appropriately be recast 
into (USD type, which describes a rotating axisym- 
metric SPCs. We conclude on the basis of these 
observations that the energy can be carried away 
to outside world by the thermal defuse blackbody 
X-ray radiation through a thin belt area S. As 
M82X-2 spins, the pulse profile results from the 
axial tilt or obliquity. Hence the X-ray beams in¬ 
tercept Earth-like lighthouse beacons. The orbital 
motion causes a modulation in the observed pulse 
frequency. The SIMBH model of M82X-2 in bi¬ 
nary system is schematically plotted in Fig. [2] No 
eclipse condition holds. The parameters of a bi¬ 


nary system is viewed in the orbital plane. 


3.3. Basic geometry: Implications on the 
pulse profile and mass scaling 


To see where all this is leading to, let us consider 
next the real issue that of the physical description 
of M82X-2. A knowledge of the dynamical mass 
measurements of the compact objects that power 
ULXs is a primarily necessary prerequisite to the 
derivation of a complete picture about the phys¬ 
ical nature of ULXs. Keeping in mind aforesaid, 
we are now in a position to derive a general pulse 
profile dependent upon the position angles, and 
give a quantitative account of a potential dynam¬ 
ical mass scaling of M82X-2 and other energetics. 
The most reliable method is to measure the mass 
function through the secondary mass and orbital 
parameters, which can be measured only if the sec¬ 
ondary donor star is optically identified. In the ab¬ 
sence of direct mass-function measurements from 
phase-resolved optical spectroscopy, we still have 
to rely on X-ray spectral and timing modeling and 
other indirect clues. To bring this point sharply 
into focus note that in case of the first ultralumi- 
nous pulsar, only the X-ray mass function is mea¬ 
sured ( Bachetti et al. 112014 1. when the optical sec¬ 
ondary is unknown and most of the orbital param¬ 
eters are yet to be measured. However, exploring 
the key physical characteristics of a SIMBH model, 
let us consider the space-fixed Cartesian coordi¬ 
nate system labeled (z,x,y), with zx as a plane-of 
sight, and the axis s of the M82X-2-fixed frame 
as the rotation axis. A schematic plot is given in 
Fig. El Here and throughout we now use following 
notational conventions. The angles 6 and (j) are 
spherical polar coordinates. The observed pulses 
are produced because of periodic variations with 
time of the projection on the plane-of sight, dzx(t)^ 
of the vector d{t) collinear to n{t) {d{t) = d ), 


where n{t) is the normal to the plane of the proto¬ 
matter disk at the moment t. The n(0) lies in the 
plane of zs. So, these pulsations are due to the 
fact that the visible surface is less at one moment 
than at another. The Fig. El explains the spherical 
triangle. That is, given a unit sphere, a spherical 
triangle on the surface of the sphere is defined by 
the great circles connecting three points u, v, and 
w on the sphere (shown at top). The lengths of 
these three sides are a = (s,n) (from u to v), the 
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Fig. 3.— A schematic plot explaining the spherical triangle solved by the law of cosines. The space-fixed 
Cartesian coordinate system is labeled (z,x,y), with zx as a plane-of sight. Axis s of the M82X-2-fixed frame 
is rotation axis. The angles 0 and </> are spherical polar coordinates. The line of nodes N is defined as the 
intersection of the equatorial and proto-matter disk planes. It is perpendicular to both the z axis and vector 
n{t), where n(t) is the normal to the plane of proto-matter disk at moment t. The n(0) lies in the plane of 
zs. The lengths of three sides of a spherical triangle (shown at top) are 0 = {z,s), a = (s,n) and P = {z,n). 
The vertex angle of opposite P is Vlt. 


axial tilt 9 = {z,s) (from u to w), and p = {z,n) 
(from V to w). The angles of the corners u and e 
opposite P equal u = e = fit. The proto-matter 
disk was shifted from the orbital direction on an¬ 
gle 6 = 9 — a towards the equator. The projection 
dzx{t) is then written 

dzx{d, 9, 4>, a, t) = a / dl{d, 9, a, t) + dl{d, 0, a, p, t), 

(14) 

where 

dz{d,9,a,t) = d{...,t) ■ el = dcosP{9,a,t), 

(15) 

and 

dxid,9,a,(j),t) = d{...,t) ■ el = 
d sinP{9,a,t) cos (j)n{9,a,4>,t). 

Here the Cz and Cx denote unit vectors along the 
axes z and x, respectively, </>„ = ^ -I- A is the 
azimuthal angle of vector n{t). The vertex an¬ 
gle opposite the side a is A. To the extent that 
all of the rotational energy of M82X-2 is dissi¬ 
pated as thermal defuse X-ray blackbody radia¬ 
tion, this may escape from the event horizon to 
outside world only through a thin belt area S. 
The radiation arisen from per area of surface is 
ctT^, where is the surface temperature, <7 is the 


Stefan-Boltzmann constant. Therefore, the pulsed 
luminosity L will be observed if and only if the pro¬ 
jection of the belt area Szx = 27rllgdzx(d, 9, </>, a, t) 
on the plane-of sight zx is not zero. So, pulsed lu¬ 
minosity reads 

L{Rg,d,Ts,9,(j),a,t) = Szx crT^ = 
2'KRgdzx{d, 9, p, a, t) aTf = (17) 

Lo{M,d,Ts)^{9,p,a,t), 

where the amplitude and phase, respectively, are 
Lo{M,d,Ts) ~ 1.05 X lO"* (ergs"^)-^ — 

Mq m 


^{9,p,a,t) = \jf— sin^ P sin^((/() + A). 

(18) 

The spherical triangle is solved by the law of 
cosines 


cos P{9, a, t) = cos 0 cos a -I- sin 0 sin a cos fit, 

\ cos a — cos 0 cos/3 

cos A{9, a, t) = - ■ n ■ o -• 

sm 9 sm p 


( 19 ) 
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Consequently, the pulsed flux can be written in 
the form 


F{Rg,d, 9, (j), a, t) = Fo{M, d) $(6>, 0, a, t). 

( 20 ) 

where, given the distance D cs 3 .6 Mpc to the 
galaxy M82 ( Bachetti et al. 2014 1. the flux am¬ 
plitude is 


Fo{M,d) 


Lo{M,d) 

FkD'^ 


6.8x10 '^®(ergs ^cm 


M d 


( 21 ) 


Thus, the theoretical model of periodic source 
M82X-2 left six free parameters: {M,d,Ts,6, </), a). 
The figure Fig. |4] reveals the diversity of the be¬ 
havior of characteristic phase ^{9,(j),a,x = ilt) 
profiles versus the time, viewed at given position 
angles {9, (j), a). The observed X-ray pulse is fur¬ 
ther determined by the complicated transfer of X- 
ray photons from the surface of M82X-2 through 
regions of external accreting plasma. At hard 
look, the position angles being the parameters of 
a model function can be evaluated via nonlinear 
regression analysis to the approximate solution of 
overdetermined systems to best fit a data set from 
observed pulsed profile of M82X-2. These missing 
ingredients are the shortcoming of present frame¬ 
work, which will be solved by iterative refinement 
elsewhere. Now according to (HZl, for maximum 
value of pulsed luminosity either at /3 = tts or 
(j) + A = TTs (s=0,l,2,...), we have Lo{M,d,Ts) = 
Z(3 - 30keV) = 4.9 x lO^^ergs-i or Fo{M,d) = 
F{3 — 30keV) ~ 3.16 x 10“^^ ergs“^cm“^. Hence 
the surface temperature scales oc with the 
black hole mass: 

M _ 4.66 X 10^^ m 

Tj 7 ’ 


That is, a cooler radiation surface implies a big¬ 
ger black hole. If we further assume that the 
persistent emission -L(0.3 — lOkeV) = 1.8 x 
104 Oergg-i- from M82X-2 is isotropic, we may 
impose a strict Eddington limit on the mass trans¬ 
fer rate that can be accepted by the black hole, 
L < LEdd — 1-3 X 10^® ergs“^. This imposes 
stringent constraint on the lower limit of black hole 
mass 

M 

— > 138.5, Rg > 408.6 km, (23) 


because for an accreting of ~ 10 per cent of the Ed¬ 
dington limit - a fairly typical accretion rate for a 
high-state black hole - this points towards a rough 
estimate of upper limit of mass M < 1385 Mq. 


3.4. More accurate determination of up¬ 
per mass limit of M82X-2 

For the knowledge of upper mass limit of M82X- 
2 , a good progress can still be made by estab¬ 
lishing a direct physical relation between masses 
of M82X-2 and M82X-1, and then one should 


fe.e. lOkaiima et al. (2006h: 

Mucciarelli et al. 

(2006ll: Dewangan et al. (12006 

1: Stobbart et al. 

(2006[l: 

Patruno et al. (200611: Portegies et al. 

(200i; 

Roberts 1 (20071: Casella et al. ( 

20081: 

Feng & Kaaret 1 (2010ll: Feng & Soria ( 

20111; 


Pasham et al. { 2014^ '). The controversy, however 


is with their mass range. As one may envisage, 
a different mass estimates of M82X-1 may yield 
different mass values of M82X-2. So, we should 
be careful in choosing the most accurate black 
hole mass measurement to-date. As the centroid 
of the persistent emission is between M82X-2 and 
M82X-1 which indicates that M82X-2 harbors an 
IMBH, likewise M82X-1, we suppose that accre¬ 
tion onto a black hole is well approximated by 
the relation Lace = rjc^M = This gives 

Ml/M = Li/L ~ 5.556. According to MTBH, 
we have M oc for both the collisionless and 
the h ydrodynamic spherical accr etions onto black 
hole ( Ter-Kazarian 20l4 2015a 1. Making use of 
these relations gives 


Ml ~ 2.357M. 


(24) 


The M82X-1 is a good candidate for hard state 
ULXs which may be one of the very few ULXs 
that change their spectral state during out¬ 
bursts, switching from a hard to a thermal 


state (jFeng fc Soria 1120111 1. The type-C low fre¬ 


quency quasi-periodic oscillations and broadband 
timing noise, detected in the two XMM-Newton 


observations in 2001 and 2004, in the central 


region of M82 (IStrohmaver fc Mushotzkv 


Dewangan et al. 


2006 t 


Mucciarelli et al. 



an d later confirmed to o riginate from the M82X- 


1 (|Feng fc Kaaret 1200711 . suggest that the ULX 
harbors a massive black hole. The mass esti¬ 


mate by (jPewangan et al. 1120061 1 is based on the 
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Fig. 4.— Characteristic phase profiles versus the time {x = nt), viewed at the following position angles 
(1) (45°,30°,0°); (2) (45°, 30°, 90°); (3) (45°, 30°, 135°); (4) (45°, 30°, 60°); (5) (45°, 90°, 90°); 
(6) (45°, 18°, 90°); (7) (45°, 60°, 90°); (8) (45°, 78.3°, 90°); (9) (90°, 90°, 180°); (10) (90°, 0.6°, 72°); 

(11) (90°, 30°, 60°); (12) (72°, 30°, 60°); (13) (60°, 90°, 60°); (14) (60°, 45°, 60°); (15) (60°, 153°, 60°); 

(16) (1.2°, 90°, 60°). 


assumption that M82X-1 also follows well estab¬ 
lished relation of the photon spectral index versus 
QPO frequency, F — vqpo^ found for the Galac¬ 
tic X-ray binaries in their high or intermediate 
states. The resulting IMBH mass for M82X-1 
is in the range of 25 — 520 Mq. However, there 
may be systematic errors in the photon indices 
measured with XMM-Newton and Rossi X-ray 
Timing Explorer (RXTE) due to contamination 
from nearby sources, as indicated by large appar¬ 
ent changes in the effecti ve absorption colum n. 
Another mass estimate by ( Casella et al. 1 [200i) is 
inferred from the correlations with the X-ray lumi¬ 
nosity and type C QPO frequency. This method 
is based on the correlation between character¬ 
istic frequencies, on the fundamental plane and 
on the variability plane of accreting black holes. 
Exploring this method, the black hole masses in¬ 
ferred from the characteristic frequencies are all 
about 1 0° — 10"*^ Mirt indicating that ULXs contai n 
IMBHs (|Casella et al. II2008I: iFeng fc Soria 1M) . 
But, these results were not without problems, no¬ 
tably pointed out by the same authors. Such mass 
estimates are based on scaling relations which use 


low-frequency characteristic timescales which have 
large intrinsic uncertainties. In particular, it was 
unclear whether these mHz oscillations are indeed 
the Type-C analogs of stellar mass black, and both 
the Type-C and the mHz oscillations are variable, 
resulting in a large dispersion in the measured 
mass of 25 — 1300 Mq. Positive identification 
of the emission states requires both timing and 
spectral information. Consequently, with simul¬ 
taneous observations exploiting the high angular 
resolution of Chandra to isolate the ULX spec¬ 
trum from diffuse emission and nearby sources 
and the large collecting area of XMM-Newton 
observ ations of M82 to obta in timing informa¬ 


tion, Feng fc Kaaret I (20 inl l found surprisingly 


that the previously known QPOs in the source 
disappeared. The light curve was no longer highly 
variable and the power spectrum was consistent 
with that of white noise. The energy spectrum also 
changed dramatically from a straight power-law to 
a disk-like spectrum. The disappearance of QPOs 
and the low noise level suggest that the source 
was not in the hard state. All results are well con¬ 
sistent with that expected for the thermal state. 
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The monitoring data from RXTE indicate that 
these Chandra and XMM-Newton observations 
were made during the source outbursts, suggesting 
that M82X-1 usually stays in the hard state and 
could transition to the thermal state during out¬ 
bursts. The spectral fitting suggests that the ULX 
contains a close to Eddington {Ldisk/LEdd ^ 2) 
rapidly spinning IMBH of 200 — 800 Mq masses. 
The thermal dominant states are all found dur¬ 
ing outbursts. Nonetheless, modeling of X-ray 
energy spectra during the thermal-dominant state 
using a fully relativistic multi-colored disk model 
has large uncertainties owing to both systematic 
and measurement errors. In addition to the large 
mass uncertainty associated with the modeling, 
the same authors also found that the energy spec¬ 
tra can be equally well-fit with a stellar-mass black 
hole accreting at a rate of roughly 160 higher 
than the Eddington limit. Also, the X-rays from 
this source are known to modulate with an or¬ 
bital periodicity of 62 days, which indicates to 
an intermediate-mass black hole with mass in the 
range of 200 — SOOO Mrr^ ( Patruno et o.l. I 120061: 
Portegies et al. l[2007l l. But, a recent study finds 
that this periodicity may instead be due to a pro¬ 
cessing accretion disk in which case a stellar-mass 
black hole will s uffice to explain the apparen t 
long periodicity ( Pasham fc Strohmaver I 1201, ih . 
Thus, the mass measurements above have large 
uncertainties. This makes black hole masses ob¬ 
tained by this method at the very least ques¬ 
tionable. In what follows, therefore, we adopt 
an alternative and less ambiguous mass deter¬ 
mination f or intermediate-mass black holes sug¬ 
gested by ( Pasham et al. I liHl which seems to 
be a more reliable determinant of the mass of 
M82X-1. These authors searched RXTE’s pro¬ 
portional counter array archival data to look for 
3 : 2 oscillation pairs in the frequency range of 
1 — 16 Hz which corresponds to a black hole 
mass range of 50 — 2000 Mq. In stellar-mass 
black holes it is known that the high frequency 
quasi-periodic oscillations that occur in a 3 : 2 
ratio (100 — 450 Hz) are stable and scale in¬ 
versely with blac k hole mass with a reasonabh 


verseiy witn blac k note mass witn a reasonablY 
small dispersion (iMcClintock fc RemillarcT 200^ 


Remillard et al fooSll 2002aibl : Strohmaver 2001al f9l. 

Pasham et al. I 1201411 report stable, twin-peak 


3.32 ±0.06 Hz (coherence, Q = centroid frequency 
(i/)/width(!/) > 27) and 5.07 ±0.06 Hz (Q > 40). 
The discovery of a stable 3 : 2 high-frequency pe¬ 
riodicity simultaneously with the low-frequency 
mHz oscillations allows for the first time to set 
the overall frequency scale of the X-ray power 
spectrum. This result not only asserts that the 
mHz quasi-periodic oscillations of M82X-1 are the 
Type-C analogs of stellar-mass black holes but also 
provides an independent and the most accurate 
black hole mass measurement to-date. Assuming 
that one can scale the stellar-mass relationship, 
they estimate the black hole mass of M82X-1 to 
be 428 ± 105 Mq. They also estimate the mass us¬ 
ing the relativistic precession model, which yields 
a value of 415 ± 63 Mq. C ombining the average 2 
10 keV X-ray luminosity Pasham fc Strohmaver 


(2013); Kaaret fc Peng ( 2007tl of the source of 
5 X lO'^^ergss”^ with the measured mass suggests 
that the source is accreting close to the Eddington 
limit with an accretion efficiency of 0.8 ± 0.2. 

Making use of the mass values 428 ± 105 Mq 
with ((24)l . we provide, therefore, the mass estimate 
for M82X-2: 

M- 138.5 - 226 Mq, i?g ~ 408.6 - 666.7km, 

(25) 

Rotation speed at surface of M82X-2 with upper 
limit mass 226 Mq as rigidly spinning IMBH con¬ 
figuration of angular velocity H equals w = ~ 

3.06 X 10®cms“^. We also have slightly improved 
the lower mass limit 32 3 Me, of M82X-1 given 


by (jPasham et al. 


20141 1 now to be 326.5 Mq. 
we obtain then 


Combining ([25l) and 


7^4 j 

2.06 X 10^3 < — < 3.34 X 10' 

m 


,33 


(26) 


(3:2 frequency ratio) X-ray quasi-periodic os¬ 
cillations from M82X-1 at the frequencies of 


Eor reasons that will become clear below, next, in 
our setting we retain the rat her concrete proposal 
of non -spinning black holes (ITer-Kazarian 12014 , 
2015bl l. i.e. the neutrino flux from spinning M82X- 
2 might as well be highly beamed along the plane 
of proto-matter disk and collimated in very small 
opening angle. For the values (1^ . this yields the 
constraints 

7.5 X 10-'^ — < - ed < 1.2 X 10-® —. (27) 

m m 

Besides, the ed is likely to be about an order of 
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magnitude ~ 10 ^ for M82X-2. Therefore, 

d ~ 61 - 100 m, Ed (4.6 - 7.5) x 10-^ 

(28) 

is a good guess for the thickness of proto-matter 
disk at the edge of event horizon. This, together 
with (IMl), give Ts ~ 7.6 x 10^ K. Thus, M82X- 
2 indeed releases ~ 99.6% of its radiative energy 
predominantly in the X-ray bandpass of 0.3 — 30 
keV. However, the studies in other wavelengths 
well give us useful information on its physical na¬ 
ture and environment. From Wien’s displacement 
law we obtain the wavelength Xmax — 0.381nm 
at which the radiation curve peaks, which corre¬ 
sponds to energy hvmax — 3.2 keV. As an immedi¬ 
ate corollary to the assumption that the emission 
arisen from accretion is isotropic, we are able to 
infer the most important ratios of the pulsed and 
persistent luminosities to the isotropic Eddington 
limit for M82X-2: 

-^-0.17-0.28, -^- 0.63- 1.03, 

^Edd ^Edd 

(29) 

respectively, where Lsdd — (1.75 — 2.85) x 
10^0 erg s“i. These properties appear consistent 
with the sub-Eddington hard state, which given 
the observed luminosities of M82X-2 suggests the 
presence of SIMBH with a dynamical mass (Ell). 
Given the angular velocity H = = 1.37s), 

we may calculate the rotational kinetic energy 
E^ot = of M82X-2, where I = is the 

moment of inertia if M82X-2 is regarded as the 
rigidly spinning canonical configuration of mass 
M and radius Rg. Hence 

~ (3.72 - 16.17) X lO'^i erg. (30) 

When all energy thermalized, radiation emerges as 
a blackbody. A significant fraction of the accretor 
M82X-2 surface radiates the accretion luminosity 
at temperature 

such that Tb ~ (3.88 — 5.37) x 10®AT. The grav¬ 
itational energy of each accreted electron-proton 
pair turned directly into heat at (shock) tempera¬ 
ture Tsh-. SkTsh = so Tsh 1.8 x lO^^K. 

Hence typical photon energies of persistent radia¬ 
tion lies between fcTf, cs 0.34 keV < hi/ < kTsh — 


180 MeV. So, M82X-2 is persistent X-ray and pos¬ 
sibly gamma-ray emitter. Also, given the mass of 
the most brightest source M82X-1 of persistent X- 
ray radiation, typical photon energies of persistent 
radiation lie in range kTbi^xi) — 0.3 keV < hv < 
kT^f^ Ed 180 MeV. 


3.5. The mass of companion star and orbit 
parameters 

Once the mass scaling of M82X-2 is accom¬ 
plished, this can potentially be used further to 
quantify the association between the M82X-2 and 
the optical secondary donor star in X-ray bi¬ 
nary. The orbital period Por is a key parameter 
for dynamical mass measurement. Erom Eig. [51 
the separation of the two masses is a, and their 
distances from the center of mass are oi and 
02. The highly circular orbit, combined with the 
mass function f{M, M 2 , i) = 2.1Mq, the lack of 
eclipses and assumption of a Roche-lobe- filling 
compani on constrain the incl ination angle to be 
i < 60° ( Bachetti et al. 2014 ). These alow to de¬ 
termine the mass M 2 of donor star: 


M2 

Mp, 


> 


48.3, for M= 138.5 M0, 
64.9, for M = 226 Mq. 


(32) 


Thus, optical companion is a typical 0/B super¬ 
giant, which evolves away from the main sequence 
in just a few million years. The binary separa¬ 
tion can be computed with Kepler’s Law a® = 
The Doppler curve of the spectrum 


G(M+M2) o 2 


of NuSTAR J095551-b6940.8 shows a Por = 2.5- 
day sinusoidal modulatio n arisen from binary or¬ 
bit ( Bachetti et al. l[2014l ). All these give the pro¬ 
jection of the orbital velocity of M82X-2 along the 
line of sight vi = sin / cs 200.9^, and hence 

oisini cs 9.97?©. The absence of eclipses implies 
cos i > ^, where R is the radius of companion 
donor star. Hence 


R > 


22.17?©, for M= 138.5 Mq, 
25.7 7?©, for M = 226M©. 


(33) 


As well as oi > 11.4 7?©, 02 > 32.9 7?© for 
M/Mq = 138.5; and ai is the same, 02 > 39.9 7?© 
for M/Mq = 226. The Roche lobe radius Rl for 
donor star is 


Rl > 


15 7?©, for M = 138.5 Mq, 
14.3 7?©, for M = 226Mq. 


( 34 ) 
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Thus, the M82X-2 and donor star constitute the 
semi-detached binary, accreting through Roche- 
lobe overflow. Donor star exceeds its Roche lobe 
{R > Rl), therefore its shape is distorted because 
of mass transfer from donor star through the inner 
Lagrange point Li to the M82X-2. The accretion 
stream is expected to be rather narrow as it flows 
through the Li point and into the Roche lobe of 
the primary. 


3.6. Spin-up rate and the torque added to 
M82X-2 

All accreting pulsars show stochastic variations 
in their spin frequencies and luminosity, includ¬ 
ing those displaying secular spin-up or spin-down 
on long time scales, blurring the conventional 
disti nction between disk-f ed and wind-fed bina¬ 
ries ( Bildsten et al. ][1993). Pulsed flux and ac¬ 
cretion torque are strongly correlated in outbursts 
of transient accreting pulsars, but uncorrelated, or 
even anticorrelated, in persistent sources. The ob¬ 
served secular spin-up rate can be accounted for 
quantitatively if one assumes the reduction of the 
torque on the rapidly spinning object. Continu¬ 
ing on our quest, below we determine the condi¬ 
tions under which pulsed source M82X-2 spins up 
and hance gains rotational energy as matter is ac¬ 
creted, i.e. we discuss the relationship between the 
properties of the exterior flow and the measured 
rate of change of angular velocity dVl jdt. Although 
measurements of spin-up or spin-down appears to 
be the most promising method for determining 
the angular momentum transport by the inflowing 
matter, which in turn, may provide information 
about pattern of material flow outside the event 
horizon of SIMBH, the extraction of this infor¬ 
mation from such measurements clearly requires 
some care. We explore the relationship between 
the torque (= 1 ) flux through the event horizon, 
spin-up rate of SIMBH and the rate of change of 
its rotational energy. The rates of change of the 
SIMBH angular velocity and of the rotational en¬ 
ergy can be related to the flux of torque across 
the event horizon boundary as foll ows. The rate 


of ch ange of the torque is given by (IGhosh et al. 
197711 

j^{m)=Ml, (35) 

where I is the moment of inertia, and I is the 
torque added to the SIMBH per unit mass of ac¬ 


creted matter. Equation (1351) gives for the rate of 
change of angular velocity 

dn _ M 

dt Lbh 


in - n^Ri 


M di \ 

Tm) 


(36) 


where L^h = m, and Rgtr is the radius of gyration 
of SIMBH. The rate of change of the rotational 
energy is 


= A f ImA . 

dt dt \2 J 


(37) 


To make further progress we recast dMl) and (EZD 
into the form 


dEj‘ot 

dt 


= M 



2 

gir 


fM dl \ 

[r'dM) 


(38) 


When > 0, which is generally the case, the 

SIMBH’s behavior can be conveniently character¬ 
ized by the dimensionless parameter 


c = 


I 


/M dl\~^ 
\TdM) 


(39) 


Thus the black hole loses rotational energy and 
spins down for ^ < 1/2, whereas it gains rotational 
energy and spins up for ^ > 1; for 1/2 < C, < 1 
the black hole spins down even though it is gain¬ 
ing rotational energy. The logarithmic derivative 
canonical configuration, i.e. spinning 
uniform-density sphere with mass M and radius 
Rn , is 


din J 
dlnM 


dlnM 


In (^^MRg^ = 3, so 


c 


I 

3Di^’ 


(40) 


where Rg^^ = I/M = For the spin up regime 

of M82X-2 when C ~ 0.073 > 1, we obtain 

^gir 

I > 2.192Rls-\ so 


J 3.7 X 10'^ km^s-^ for M=138.5Mq 
^ \ 9.7 X lO^km^s-i, for M = 226Mq. 

(41) 

The time derivative of the angular velocity (1361) 
gives 


= ^ 3 
5P 


P m\ 

Da ?j ■ 


(42) 


Combing M ~ 6.35 x 10 ^Mq yr and a linear 
spin-up p of the NuSTAR J095551-I-6940.8 pulsar. 
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from (1421) we obtain the torques added to M82X- 
2 per unit mass of accreted matter, which satisfy 
the spin-up condition (IdTll of C > 1: 


I ~ 


1.1 X lO^km^s-i, for M = 138.5M©, 
9.8 X lO^km^s-i, for M = 226Mq. 

(43) 

This is not a final report on a closed subject, and 
needless to say once again that astronomers are 
not yet at the end of producing a coherent picture 
of the ULX Universe, so there is plenty to argue 
about. We believe that there is still very much 
to be gained by further study of the issues that 
we raised in this paper. There are deep concep¬ 
tual and technical problems involved, and these 
provide scope for the arguments discussed, which 
are carefully presented in both mathematical and 
physical terms. It should be emphasized that the 
key to our construction procedure of both MTBH 
in general, and suggested model of mighty ultra- 
luminous X-ray pulsations from periodic source 
M82X-2 in particular, is widely based on the 
premises of our experience of accretion physics. 
Therefore, what we have presented here has all 
the vices and virtues of the classical scenario of 
runaway core collapse which has always been a 
matter of uncertainties and controversies. On the 
other sight, astronomers doubted the existence of 
the superdense proto-matter sources, away from 
the accretion physics. Models of this type have 
a long tradition of precursors th at dates back to 
the pioneering sem inal papers by Ambartsiimian I 
( 1958. 1961. 1965 1. Since then the very possibility 
of existence of such sources was rejected, and the 
idea of the superdense proto-matter was dismissed 
as a fancy of eccentric astronomer. However, we 
caution that our entire constructions will be valid 
as well in the case if some hitherto unknown yet 
mechanism in Nature will in somehow or other way 
produce the superdense proto-matter, away from 
the accretion physics. 


4. Concluding remarks 

Our analysis of periodic source M82X-2 goes 
quite on the contrary to the conventional wisdom. 
The pulsed luminosity of M82X-2 is the most ex¬ 
treme violation of the Eddington limit and could 
be reconciled with that in the model of common 
pulsar only by very arbitrary assumptions on geo¬ 
metric beaming of accretion flow on neutron star. 


Instead of making such assumptions, our approach 
finds it preferable to return in part to the ideas of 
MTBH framework to circumvent the alluded ob¬ 
stacles without the need for significant breaking of 
Eddington limit. The new conceptual element of 
the implications of MTBH in tackling this problem 
is noteworthy. This explores the most important 
processes of spontaneous breaking of gravitation 
gauge symmetry and rearrangement of vacuum 
state at huge energies, making room for growth 
and merging behavior of black holes. Putting 
apart the discussion of inherent problems of the 
mass scaling of the black holes in ULXs, we have 
focused on black hole. We assume the M82X-2 is 
being SIMBH, resided in the final stage of growth. 
If this were the case, eventually, a thin co-spinning 
proto-matter disk driven by accretion would be 
formed around the spinning proto-matter core. It 
was tilted from the plane of accretion on a defi¬ 
nite angle S towards the equator and has reached 
out the edge of the event horizon, where a met¬ 
ric singularity inevitably disappears. The energy 
is carried away then from event horizon through 
a thin belt area to outside world by both the 
ultra-high energy neutrinos produced in the su¬ 
perdense proto-matter medium, and the thermal 
defuse blackbody radiation released from the outer 
surface layers of ordinary matter of spinning SPC 
and co-spinning proto-matter disk. All of the ro¬ 
tational energy of SIMBH is dissipated as ther¬ 
mal defuse blackbody X-ray radiation to outside 
world. As M82X-2 spins, we see pulses because 
of the axial tilt or obliquity. That is, the M82X- 
2 is the emitter of both persistent and pulsating 
X-rays. We derive the general prohles of pulsed lu¬ 
minosity and X-ray flux of M82X-2. Hence M82X- 
2 indeed releases 99.6% of its radiative energy 
predominantly in the X-ray bandpass of 0.3 — 30 
keV. The resulting theoretical model necessarily 
includes a number of poorly known parameters. 
We give a quantitative account of all the energet¬ 
ics, a dynamical mass scaling and orbital param¬ 
eters of the semi-detached X-ray binary contain¬ 
ing primary M82X-2 and the secondary massive 
0/B-type donor star, accreting through Roche- 
lobe overflow. The position angles can be evalu¬ 
ated from rigorous comparison with the behavior 
of observed pulsed light curve of M82X-2, which 
will be discussed in a future publication. 
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